Constraints of Relic Gravitational Waves by Pulsar Timing Array: Forecasts for the 
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Measurement of the pulsar timing residuals provides a direct way to detect relic gravitational 
waves at the frequency / ~ 1/yr. In this paper, we investigate the constraints on the inflationary 
parameters, the tensor-to-scalar ratio r and the tensor spectral index nt, by the current and future 
Pulsar Timing Arrays (PTAs) . We find that Five-hundred-meter Aperture Spherical radio Telescope 
(FAST) in China and the planned Square Kilometer Array (SKA) projects have the fairly strong 
abilities to test the phantom-like inflationary models. If r = 0.1, FAST could give the constraint 
on the spectral index nt < 0.56, and SKA gives nt < 0.32. While an observation with the total 
time T = 20yr, the pulsar noise level a w = 30ns and the monitored pulsar number n = 200, could 
even constrain n t < 0.07. These are much tighter than those inferred from the current results of 
Parkers Pulsar Timing Array (PPTA) , European Pulsar Timing Array (EPTA) and North American 
Nanohertz Observatory for Gravitational waves (NANOGrav). Especially, by studying the effects 
of various observational factors on the sensitivities of PTAs, we found that compared with a w and 
n, the total observation time T has the most significant effect. 

PACS numbers: 04.30.-w, 04.80.Nn, 98.80.Cq 



I. INTRODUCTION 

In a whole range of scenarios of the early Universe, in- 
cluding the well-studied inflationary models, a stochas- 
tic background of relic (primordial) gravitational waves 
(RGWs) was produced due to the supcradiabatic ampli- 
fication of zeropoint quantum fluctuations of the grav- 
itational field [1H3. Their detection maybe provide the 
unique way to study the birth of the Universe, the expan- 
sion history of Universe before the recombination stage, 
and test the applicability of general relativity and quan- 
tum mechanics in the extremely high-energy scale [J] . 

Since RGWs have a wide range spreading spectra, from 
10 _18 Hz to 10 10 Hz, one can detect or constrain them at 
different frequencies. The temperature and polarization 
anisotropies of Cosmic Microwave Background (CMB) 
radiation provide the way to constrain RGWs at very low 
frequencies, / < 10~ 15 Hz. Nowadays, combining with 
other cosmological observations, the nine- year WMAP 
data place the constraint on the tensor-to-scalar ratio 
r < 0.13 [![. While the new Planck data give the the 
tightest constraint r < 0.11 6], which is equivalent to 
the constraint of the amplitude of RGWs at lowest fre- 
quency / ~ 10~ 17 Hz. In the near future, this bound 
will be greatly improved by the forthcoming polarization 
observations of Planck satellite, several ground-based 
and balloon-borne experiments (BICEP, QUIET, PO- 
LARBEAR, QULIOTE, ACTPOL, SPTPOL, QUBIC, 
EBEX, PIPER, SPIDER ct al), and the planned fourth- 
generation CMB missions (CMBPol, LitcBird, COrE. ct 
al.). 

Among all the direct observations, LIGO S5 has ex- 
perimentally obtained so far the most stringent bound 



n gw < 6.9 x 10~ 6 around / - 100Hz @, Q. It is ex- 
pected that AdvLIGO, AdvVIRGO, KAGRA, ET and 
eLISA will also deeply improve it in the near future. In 
particular, the Planned BBO, DECIGO and ASTROD 
projects may directly detect the signal of RGWs in the 
far future. In addition, there are two bounds on the in- 
tegration J Q gw (f)dlnf <, 1.5 x 10 -5 , obtained by the 
Big Bang nucleosynthesis (BBN) observation Q and the 
CMB observation [lo[. 

By analyzing of pulsar pulse time-of-arrival (TOA) 
data, people find the millisecond pulsars are very sta- 
ble clocks. The measurement of their timing residuals 
provides a direct way to detect GW background in the 
frequency range / € (10~ 9 , 10~ 7 )Hz [IlHlJ]. In addi- 
tion to the GWs generated by the coalescence of massive 
black hole binary systems [3] and cosmic strings fl5j |. 
RGWs are another kind of most important GW sources 
in this frequency range. Recently, PPTA, EPTA and 
NANOGrav teams have reported their observational re- 
sults on the stochastic background of GWs. In [l6j|, by 
considering these results, wc have dctailcdly investigated 
the constraints on the Hubble parameter during inflation 
in the most general scenario for the early Universe. In 
this paper, wc shall extend them to constrain the tensor- 
to-scalar ratio r and the tensor spectral index n t . 

In addition, as the main goal of this paper, we will 
discuss the potential constraint (or detection) of RGWs 
by the future PTA observations. In our discussion, FAST 
and SKA will be treated as two typical projects, and 
mainly focused on in the studies. The dependence of the 
RGWs constraints on the total observation time T, the 
number of monitored pulsars n and the magnitude of the 
pulsar timing noise o~ w will be discussed. 
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This paper is constructed as follows, In Sec. 2, wc 
briefly review the model to describe the RGWs, and re- 
late the energy density of GWs il gw , the tensor-to-scalar 
ratio r and the tensor spectral index n t to the character- 
istic strain spectrum h c (f), which is widely used in the 
PTA analysis. In Sec. 3, we describe the sensitivities 
of the current and future experiments, and discuss the 
dependence on various observational parameters. Sec. 4 
summaries the main results of this paper. 



II. RELIC GRAVITATIONAL WAVES 

Incorporating the perturbation to the spatially flat 
Fricdmann-Robertson- Walker spacetime, the metric is 



ds 2 = a 2 {rf) \drj 2 — (Sij + hijdx l dx-')\ , 



(1) 



where a is the scale factor of the universe, and 7/ is 
the conformal time, which relates to the cosmic time 
by adrj = dt. The perturbation of spacetime fty is a 
3x3 symmetric matrix. The gravitational-wave field is 
the tensorial portion of hij , which is transverse-traceless 



dih? 3 = 0, S^hi 



0. 



RGWs satisfy the linearized evolution equation [l|: 
-gd^hij) = —l&ivG-Kij. 



(2) 



The anisotropic portion 7Ty- is the source term, which 
can be given by the relativistic free-streaming gas fl7| | . 
However, it has been deeply discussed that the rela- 
tivistic free-streaming gas, such as the decoupled neu- 
trino, can only affect the RGWs at the frequency range 
/ € (1CT 16 , l(T 10 )Hz, which could be detected by the 
future CMB observations So, it cannot obviously in- 
fluence the RGWs at the frequency / € (1CT 9 , 10~ 7 )Hz. 
For this reason, in this paper we shall ignore the contri- 
bution of the external sources. So the evolution of RGWs 
only depends on the scale factor and its time derivative. 

It is convenient to Fourier transform the equation as 
follows: 



d 3 k 
(2tt) 3 / 2 



E 

r=+,X 



(s) As) ik-S 



+ C.C. 



(3) 

where c.c. stands for the complex conjugate term. The 
polarization tensors are symmetry, transverse-traceless 

0, and satisfy the conditions 



k l e[ s ){k) = 0, S^efMk) 



eWfyeWtf) = 26 ss , and e\ s /(-k) = e\ s /(k). Since the 
RGWs we will consider are isotropy, and each polariza- 
tion state is the same, we have denoted hrf\ri) by hk(n), 

where k = \k\ is the wavenumber of the GWs, which re- 
lates to the frequency by k = 2nf. (The present scale 
factor is set oq = 1). So Eq.([2|) can be rewritten as 



h k " + 2-h k ' + k 2 h k =0, 
a 



(4) 



where the prime indicates a conformal time derivative 
d/dr]. For a given wavenumber k and a given time ry, we 
can define the transfer function tf as 



t f (r),k) = h k (r))/h k {r)i), 



(5) 



where rji is the initial time. This transfer function can 
be obtained by solving the evolution equation 

The strength of the GWs is characterized by the 
GW energy spectrum, f2 gw = p gw /po, where p gw = 



;{hijh % i), the critical density is po 



3Hx 



is the current Hubble constant. Using Equations in fl5|) 
and §5$ , the energy density of GWs can be written as |19| 



Pgw 



dkPt{k)i 2 ( m ,k) 



32ttG 



(6) 



where Pt(k) = ^pr\hk(r]i)\ 2 is the so-called primordial 
power spectrum of RGWs. Thus, wc derive that the cur- 
rent energy density of RGWs, 



gw 



/ 



fl gw (k)d\nk, and r2g W (fc) 



Pt(k); 



12H§ 



tf(Vo,k), 



(7) 

where the dot indicates a cosmic time derivative dj dt. 

Now, let us discuss the terms Pt(k) and tf(r/o, k) sepa- 
rately. The primordial power spectrum of RGWs is usu- 
ally assumed to be power-law as follows: 



(8) 



P t (k) = A t (K) [ - 



This is a generic prediction of a wide range of scenar- 
ios of the early Universe, including the inflation models. 

At(k*) = ^—r directly relates to the value of Hubble pa- 
rameter H at time when wavelengths corresponding to 
the wavenumber fc* crossed the horizon fll. Il6l l2(ij. In 
observations, wc always define the tensor-to-scalar ratio 
r, and write the amplitude of RGWs as At{k*) = A s r, 
where A s is the amplitude of primordial density pertur- 
bation at k = k*. nt is the spectral index of RGWs, 
which relates to the effective equation-of-state w of the 
cosmic "matter" in the inflationary stage by the relation, 



n t 



l + 3u> 



(9) 



If the inflation is an exact de Sitter expansion stage 
with w = — 1, we have the scale- invariant spectrum with 
n t = 0. For the canonical scalar-field inflationary mod- 
els, wc have w > — 1, which predicts the red spectrum 
of RGWs with n t < @. However, for the phantom 
inflationary models [2l[, one has w < — 1 and n t > 0. So 
the determination of n t can distinguish different kinds of 
inflationary scenarios. 

Now, let us turn to the transfer function tt, defined 
in ©, which describes the evolution of GWs in the ex- 
panding Universe. From Eq.@, wc find that this transfer 
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function can be directly derived, so long as the scale fac- 
tor as a function of time is given [22l42q . In this paper, 
we shall use the following analytical approximation for 
this transfer function. It has been known that, during 
the expansion of the Universe, the mode function hk(rj) 
of the GWs behaves differently in two regions [22| . When 
waves are far outside the horizon, i.e. k <C aH , the ampli- 
tude of hk keeps constant, and when inside the horizon, 
i.e. k 3> aH, the amplitude is damping with the expan- 
sion of Universe, i.e., oc l/a(rj). In the standard hot 
big-bang cosmological model, we assume that the infla- 
tionary stage is followed by a radiation dominant stage, 
and then the matter dominant stage and the A domi- 
nant stage. In this scenario, by numerically integrating 
Eq-O, one finds that the damping function tf can be 
approximately described by the following form [27H30I ] 



-3j2(fc?7o)fi 
ki] 




1.36(— ) 



2.50(A)2 ; 
(10) 

where k eq = 0.073S! m /i 2 Mpc is the wavenumber cor- 
responding the Hubble radius at the time that mat- 
ter and radiation have equal energy density, and rjo = 
1.41 x 10 4 Mpc is the present conformal time. The factor 
f2 m encodes the damping effect due to the recent accel- 
erating expansion of the Universe [HI, H3, H3|. In this 
damping factor, we have ignored the small effects of neu- 
trino free-streamin g |l7l| and various phase transitions in 
the early Universe [26j. In this paper, we shall focus on 
the wavenumber k 3> k eq . In this range, we have the 
current density of RGWs as follows, 



15 n 2 m A s r 

16 Hl^kl q 



(11) 



which clearly presents the dependence of the RGWs on 
various cosmological parameters. 

In the PTA analysis, people always describe the GW 
background by the characteristic strain spectrum h c (f) 
[3l| . For most models of interest, it can be written as a 
power-law dependence on frequency /: 



h c (f) = A 



f 



(12) 



The characteristic strains relate to one-side power spec- 
trum P(f) and the energy density of GWs f2 gw (/) as 

h 2 (f) 2tt 2 



I2n 2 f 3 ' 



3HZ 



fK(f). (13) 



Comparing the Equations in (|11[) and (|13|). we find 
that 



«=y-l. 



and 



.4 
yr 



1 45 n? n A s 

32^ rfitf 



yr 



(14) 



(15) 



Considering the cosmological parameters based on the 
current Planck observations \m h = 0.6711, O m = 0.3175, 



Ov = 0.6825,z 



3402, A s 



2.495 x 10" y at k* 



0.002MPC" 1 [32|, we obtain that 



A = 0.88V^ x 10 5nt " 
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and 



n gw (/) = 1.09r x 10 



10nt-15 



(//yr 



-l\n t 



(16) 



(17) 



These relations will be used for the following discussion. 
Both Equations in (|16JI and (|17JI show that the amplitude 
of RGWs at / ~ 1 /yr strongly depends on the spectral 
index n t . For the cases with the scale- invariant and red 
spectrum, one always has A <5 10~ 18 and f2 gw <, 10 -15 . 
However, for the cases with blue spectrum, i.e. n t > 0, 
the values of A and i7 gw can be dramatically large. For 
example, in the models suggested by Grishchuk [33| , the 
blue spectrum with a G [—0.8, — 1] was expected, which 
corresponds to n t £ [0, 0.4], the amplitude of RGWs 



could be A ~ 10~ 17 and fi„ 



io- 



III. PULSAR TIMING ARRAY AND THE 
DETECTION OF RELIC GRAVITATIONAL 
WAVES 

A. Current constraints 

In 2006, Jenet et al. have analyzed the PPTA data and 
archival Arecibo data for several millisecond pulsars. By 
focusing on the GWs at the frequency / = 1/yr, the 
authors obtained the 2a upper limit on A as a function 
of the spectral slope a, which is presented in the left 
panel of FigfT] (black solid line) [34| . Recently, this upper 
limit has been updated by EPTA and NANOGrav teams 
[35l [3r| . It is interesting that in (34) , the authors have 
also investigated the possible upper limit (or a defini- 
tive detection) of stochastic background of GWs by using 
the potential completed PPTA data-sets (20 pulsars with 
an rms timing residual of 100ns over 5 years, which is 
also expected the case for future EPTA and NANOGrav 
projects). We have also plotted the current EPTA upper 
limit (blue dashed line), current NANOGrav upper limit 
(green dash-dotted line) and the potential PPTA upper 
limit of parameter A in the left panel of FigfT] (red dotted 
line). 

By using the relations in Eqs. (JUJ) and ([TB]) . we ob- 
tain the constraints on the parameters r and n t , which 
are presented in the right panel of FiglU Note that the 
regions above the lines are excluded by the corresponding 
PTA observations. We notice that the current Planck ob- 

which 
3j. So, 



servations give the tightest constraint r < 0.11 [6 
is nearly independent of the spectral index nt [3 
combining with Planck constraint on r, this figure shows 
the current allowed region in the r-n t plane. For example, 
if r = 0.1 is determined by the forthcoming CMB obser- 
vations, current PPTA gives the constraint n t < 0.94, 
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FIG. 1: Left panel: The upper limit of A as a function of the 
spectral slope a. Right panel: The upper limit of r for any 
given nt , where the shaded region is excluded by the current 
Planck observations. In each panel, the black solid line (i.e. 
LI) is for current PPTA 2a result [3J], the blue dashed line 
(i.e. L2) is for current EPTA 2a result [35|], the green dash- 
dotted line (i.e. L3) is for current NANOGrav 2a result [3(J 
and the red dotted line (i.e. L4) is for future PPTA 2a result 
0. 



NANOGrav gives n t < 0.90 and EPTA gives n t < 0.88 
at 2cr confident level. Meanwhile, the future PPTA will 
follow the constraint of n t < 0.67. These are listed in 
Table |U Although quite loose, these constraints would 
be helpful to exclude some inflationary models with very 
blue GW spectrum. 

From Fig Q] and Eq. (flTf . we can also obtain the con- 
straints on the energy density of RGWs fi gw (/). For 
instance, if n t = and / = 1/yr are fixed, the upper lim- 
its for fi gw (/) are listed in Table [TTl which are consistent 
with the results in [34I - I36T ] . 



For the isotropic GW background, the correlation be- 
tween the GW-induced signals are [TU, [3^, |4(j 

(44) = a 2 g H ij7kk ,, (18) 

where a g is the root mean square (RMS) of the timing 
residuals induced by GW background, which relates to 
the one-side power spectrum P(f) by a 2 g = jj^ P(f)df. 
The highest and lowest frequency of GWs are given by 
fh = and /; = ~ . Hij is the so-called Hellings-Downs 
function, which is given by Hij = lamina;— f + ^(l+8(x)), 
where x = 1 c ° b< - 6> ) and 6 is the angle distance between 
i-th and j-th pulsar, jkk' is the temporal correlation 
coefficient between the fc-th and fc'-th sampling. 

The noise term n\ includes the effects of all non-GW 
sources for the i-th pulsar. It is assumed that all noise 
sources have a fla t sp ectrum, which is consistent with 
most observations [39[ . In order to simplify the problem, 
in this paper, we assume all monitored pulsars have the 
same noise level, i.e. 

{nlK,) = alSijSkk' ■ (19) 

There are several methods to extract the GW signals 
from the observable R\ [H, |4l|, [42[. In this paper we 
follow the method suggested by Jenet et al. in 2005 j39| . 
In particular, we shall present the details of the calcula- 
tion, which are quite helpful to understand the method, 
but have been neglected in the original paper |39j ] . In ad- 
dition, some sub-dominant terms, which were neglected 
in [391, will also be presented in the finial formulas. We 
calculate the correlation coefficient between the observed 
timing residuals of each pair of observed pulsars: 

fc=l 

It is easy to get the expected values of and Cy , 

{c ij ) = o*H ij , (21) 



B. Detecting GW background by Pulsar Timing 
Array 

In the following discussion, we shall study the poten- 
tial constraints on the RGWs by the future PTA observa- 
tions, where we will focus on the Chinese FAST project 
and the planned SKA project. 

The fluctuations of the pulsar TOAs caused by the 
stochastic GW background are random. However, for 
different pulsars, these fluctuations have the correlations. 
Let us assume the observations of n 3> 1 pulsars at times 
io, t\, i Tn _i with the time interval At. The total 
observation time is T = mAt. We denote the timing 
residual of i-th pulsar at time t k as R l k , which includes 
the contribution from both GWs s\ and the noises n\, 

Le - R k = s k+ n k- 



J V m mcr^ ma^ I 

where x = Sfefe'7fcfc'/ m : an d (•) denotes the ensemble 
average. 

The comparison between Cjj and the Hellings-Downs 
function is carried out by defining the GW detection sig- 
nificance S as follows, 

v/E i - ) (««-q 2 E,- J (ff. J -«) 2 

where N = n(n — l)/2 is the number of independent 
pulsar pairs. The summation sums over all in- 

dependent pulsar pairs, i.e. ^ . = Y^h=i Sj=r The 
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TABLE I: The 2a upper limit of the spectral index n t inferred from various pulsar timing observations. 





Current PPTA 


Current EPTA 


Current NANOGrav 


Future PPTA 


FAST 


SKA 


Optimal Case 


r = 0.1 


0.94 


0.88 


0.90 


0.67 


0.56 


0.32 


0.07 


r = 0.01 


> 1 


0.99 


> 1 


0.78 


0.67 


0.44 


0.18 


r = 0.001 


> 1 


> 1 


> 1 


0.89 


0.77 


0.55 


0.31 



TABLE II: The 2a upper limit of the energy density log 10 [fl gw (/ = yr x )] inferred from various pulsar timing observations, 
where we have set n t = 0, i.e. a = — 1. 



Current PPTA 


Current EPTA 


Current NANOGrav 


Future PPTA 


FAST 


SKA 


Optimal Case 


-7.36 


-7.79 


-7.63 


-9.84 


-10.63 


-12.99 


-15.30 



quantities c and H are defined as 



i—3 



(24) 



To evaluate the quality of the detector, we need the ex- 



pected value (S), which is (S) ~ VJVcTgSjy/S c , where 
S| = ^X>«-iT) 2 > S? = l^( CjJ -(^)) 2 . (25) 



By using Eqs. ([21]) and (|22|). we get the well-known re- 
sult, 



(S) ~ VN 



X (l + iP) + ^ + ^ 
mY? H 



-1/2 



(26) 



In Jcnet et al. (2005), this formula was obtained 
by another way, which is easier to extend to the re- 
sults after low-pass filtering and whitening. It is con- 
venient to define the expected discrete power spectrum 
of R\ for the i-th pulsar P^(A,i), which includes both 
a GW component and a white noise component, i.e. 



P d (A,i) =P g (A)- 



Note that A > is the discrete 



frequency bin number corresponding to frequency A/T. 
Since we have assumed that a w has the same value for ev- 
ery pulsar, the spectrum Pd(A,i) becomes independent 
of i, so we denote it as Pj(A) in the following discussion. 
For the GW with the characteristic strain spectrum h c (f) 
in Eq. (fT2")) . one has the discrete GW- induced spectrum as 
follows, 



P (A) - ^mi^l rn(A) 
g[ >~ (2n) 2 (2-2a) [ j ' 



(27) 



where m(A = 1) = /3 2a " 2 - 1.5 2 "- 2 , and to (A > 1) = 
(A - 0.5) 2q ~ 2 - (A + 0.5) 2q " 2 . & ~ 1 is the lowest fre- 
quency used to calculate the correlation function Cy. Ac- 



tion of S r , we find that 



H 



Vp](A) + ^Vp 2 (A) 



A 



A 



^E 2 f + (l + ^ ) Vp 2(A)+(J 4 ( 



2a 



A 



")• 



and the quantity \ is calculated by ^ = P 2 (A), 
which can be gotten for any given GW background. By 
using the relation (S) ~ -y/iVcr 2 !]^ /E c , we can naturally 
obtain the result in Eq. (|2l)|) . 

In order to enhance the detection significance, the low- 
pass filtering and whitening techniques can be applied 
[39| . In this way, we can correlate only that part of signal 
which has a high signal-to-noise ratio and give each time 
series a flat spectrum to optimize the measurement of 
the correlation function. In practice, we define the new 
discrete power spectrum Pz(A) and P g (A) as follows, 



Pd(A) 



Pd(A) a\ 
P d (A) to' 



f g (A) °\ 
P d (A) m' 



(28) 



where a\ = Pd(A). In this definition, the total RMS 
fluctuation induced by GW becomes er 2 = X)a=i x 
where the summation is carried out only over the fre- 
quency bins in which the GW signal dominates the noise, 
and A max is the number of the highest frequency bin. So 
the variance S r becomes 



E 2 = ^S 2 f +^P 2 (A)+ J ff2^ P 2 (A) 



A 



A=l 



A 

° \ 



'Pi(A)' 



-4 V 2 , u d 
TO 

and the expected value of S becomes 



(29) 
,(30) 



rv A max P g (A) x 2 ^ 2 
VZ^A=1 P d (A)J H 



-1/2 



(31) 



cording to the Wicncr-Khinchin theorem and the dcfini- This formula will be used in the following subsection. 
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C. Forecasts for FAST and SKA projects 



FAST is a Chinese megascience project to build the 
largest single dish radio telescope in the world. Funding 
for FAST has been approved in 2007, and its first light is 
expected to be in 2016 [43|. It includes multibeam and 
multiband, covering a frequency range of 70MHz— 3GHz. 
The relatively low latitude (~ 26°N) of the site enables 
the observation of more southern galactic pulsars. The 
zen ith angle of FAST is about 40°, which corresponds to 
H 2 = 0.024 and Sjj = 0.155, if assuming the monitored 
millisecond pulsars evenly distribute in the observed re- 
gion. One of the scientific goals of FAST is to discover 
~ 400 new millisecond pulsars. FAST is capable of pro- 
viding the most precise observations of pulsar timing sig- 
nals, therefore, may largely increase the sensitivity of the 
spectrum window for detection of GWs. 

The noise level of the millisecond pulsars are expected 
to be a w = 30ns, after collecting the timing data for the 
total time T = 5yr [43|. As a conservative evaluate, sim- 
ilar to PPTA, we assume FAST will monitor 20 pulsars 
for the detection of GWs. Thus, by using Eq. (|3Tj) . we can 
calculate the detection significance S for any given RGW 
models, which are illustrated in FigJ5J In this figure, we 
have considered three typical models with r = 0.1, 0.01 
and 0.001. These models are predicted by the general 
inflationary models, and could be well detected by the 
future CMB observations 44|. As anticipated, if n t < 0, 
we always have 5 « 1, i.e. the detection is impossible 
for the red spectrum of RGWs. However, if the RGWs 
have the blue spectrum, the detection is possible. For 
instance, for the model with r = 0.1 and n t = 0.56 or for 
that with r = 0.01 and n t = 0.67, FAST can detect the 
signal of RGWs at 2a level. In FigEl we set (S) = 2, and 
plot the value of r for any spectral index n t . Compar- 
ing those in the right panel of FigfTJ we find that FAST 
is much more sensitive than current and future PPTA 
and/or EPTA. 

As another potential observation, we consider SKA 
project, which is a proposed major internationally- 
funded radio telescope, and is expected to be completed 
in the next decade [45j . SKA will consist of many anten- 
nas, constituting an effective collecting area of about one 
square kilometer. We expect that SKA will survey the 
full sky. If assuming the monitored millisecond pulsars 
are evenly distributed, we have H 2 = Y? H = 1/48, which 
are slightly different from those of FAST. Following [4r| . 
we assume SKA will select 100 pulsars and spend the 
total time T = lOyr for the GW detection, and the av- 
erage noise level of these pulsars are about a w — 50ns, 
which is 2 times lower than those the finial PPTA, EPTA 
or NANOGrav. In Fig|2j we consider the typical infla- 
tionary models with r = 0.1, 0.01 and 0.001, and plot 
the values of S for any n t . Again, we find the detection 
is possible, only if n t > 0, i.e. the blue GW spectrum. 
Compared with the results of FAST, the detection sig- 
nificance are much higher, due to the longer observation 
time T and the larger pulsar number n. These are also 
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FIG. 2: The detection significance of RGWs for FAST (dark 
lines, black online) and SKA (gery lines, red online) projects. 
For FAST, we have assumed T = 5yr, a w = 30ns, and n = 20, 
and for SKA we have assumed T — lOyr, a w = 50ns, and 
n = 100. For both cases, solid lines are for the models with 
r = 0.1, dashed lines are for r = 0.01, and dotted lines are for 
r = 0.001. 



clearly shown in Fig|3]and Table H In Table HI we have 
listed the detection limits of the energy density fi gw for 
the FAST and SKA projects, where we also find that 
SKA is more sensitive than FAST. 

From the formula in Eq. (|31[) . we know that the de- 
tection significance of PTA projects mainly depends on 
three factors: the total observation time T, the number 
of the monitored millisecond pulsars n and the noise level 
of the pulsar a w [47[ . Now, let us discuss the dependence 
of sensitivity on these factors separately. First, we fix 
(T w = 50ns and n — 100, and investigate the effect of ob- 
servation time T. To do it, we consider three cases with 
T = 5yr, lOyr and 20yr. Setting the detection signifi- 
cance (S) = 2, we plot the constraints of the inflationary 
models in the r-n t plane in Figf?l where we find the effect 
of total time T is very significant. For example, for the 
model with r = 0.1, the 5yr observations give the con- 
straint n t < 0.48, which can be improved to n t < 0.17 for 
the 20yr observations. For comparison, in this figure, we 
have also consider the optimal case, where a w = 30ns, 
n = 200 and T = 20yr are assumed. We find that, in 
this optimal case, the constraint of spectral index is only 
slightly improved to n t < 0.07, although the noise level 
and pulsar numbers are greatly improved. 

This effect can be understood by the following analysis. 
As well known, the contributions of GWs on the pulsar 
timing residuals mainly come from those at the lowest 
frequency range, i.e. / ~ /;. So the detection significance 
S sensitively depends on the // value. At the same time, 
we know that /; = 1/T. So the larger total observation 
time T corresponds to the smaller // value, which means 
that more low-frequency GWs can contribute the timing 
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FIG. 3: The upper limits of r and nt based on the potential 
FAST observations (black line, i.e. LI) and SKA observations 
(red line, i.e. L2). Note that the shaded region is excluded 
by current Planck observations. The dashed blue line (i.e. 
Lb) is the current tightest constraint coming from the EPTA 
observations, which is identical to that in the right panel of 
Fig. [1] The solid blue line (i.e. La) indicates the result in 
the optimal case considered in this paper, where T = 20yr, 
a w = 30ns, and n = 200 are assumed. 

residuals of pulsars. This explains why the observation 
time T is the most important factor for the sensitivity of 
PTA. 

Second, we study the effect of noise level of pulsars 
a w . Decreasing a w is equivalent to increasing the A max 
value. So a smaller a w corresponds to the case where 
more high-frequency GWs have the contributions to the 
pulsar timing residuals. However, we know that, com- 
pared with the low-frequency GWs, the high-frequency 
ones are much less important for the timing residuals. 
The results are shown in FigJSl where three cases with 
a w = 100ns, 50ns and 30ns are considered. Although as 
anticipated, lower a w corresponds to the higher sensitiv- 
ity of PTA, the effect of a w is less significant than that 
of observation time T. 

Third, the pulsar number n affects the value of S only 
by the factor \fN in Eq. ff3T|) ■ which follows that (S) cx n 
for n ^ 1. This effect is illustrated in Fig El where three 
cases with n = 50, 100 and 200 are considered. We find 
that, compared with the total observation time T and 
the pulsar noise level <j w , the pulsar number n has the 
relatively smaller influence on the detection significance 
S. 



IV. CONCLUSIONS 

Generation of GW background in the early inflation- 
ary stage is a necessity dictated by general relativity and 
quantum mechanics. The wide range spreading spectra 



FIG. 4: The upper limits of r and nt depend on the total 
observation time T. The solid black line (i.e. L2) is for the 
case with T = 5yr, dashed black line (i.e. LI) is for T = lOyr, 
and dotted black line (i.e. L3) is for T = 20yr. In all cases, 
a m = 50ns and n — 100 are assumed. The solid blue line (i.e. 
La) and dashed blue line (i.e. Lb) are identical to those in 
Fig(3] The dashed black line (i.e. LI) is identical to that for 
SKA. 

of RGWs make the possible detection at different fre- 
quency ranges by various methods. The timing studies of 
the millisecond pulsars provide a unique way to constrain 
it in the middle frequency range / G (10 -9 , 10~ 7 )Hz. 

Recently, PPTA, EPTA and NANOGrav teams have 
reported their observational results on GW background 
at / ~ 1/yx. In this paper, we infer from these bounds 
the constraint of inflation in r-nt plane. Although 
quite loose, these constraints arc helpful to exclude some 
phantom-like inflationary models. 

As the main goal of this paper, we have forecasted 
the future pulsar timing observations and the potential 
constraints on inflationary parameters r and n t , by fo- 
cusing on the FAST and SKA projects. Wc found that, 
if r = 0.1, FAST could give the constraint on the spec- 
tral index n t < 0.56, and SKA gives n t < 0.32. While 
an observation with the total time T = 20yr, the pulsar 
noise level a w = 30ns and the monitored pulsar number 
n = 200, could even constrain n t < 0.07, which can ex- 
clude or test most phantom-like inflationary models with 
this tensor-to-scalar ratio. In this paper, we have also 
studied the effects of T, a w and n on the sensitivity of 
PTA, and found that the total observation time T has 
the most important influence. So increasing the obser- 
vation time can significantly improve the sensitivities of 
the future PTAs. 
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FIG. 5: The upper limits of r and nt depend on the noise 
level a w . The solid black line (i.e. L2) is for the case with 
a w = 100ns, dashed black line (i.e. LI) is for a w = 50ns, 
and dotted black line (i.e. L3) is for a w = 30ns. In all cases, 
T — lOyr and n — 100 are assumed. The solid blue line (i.e. 
La) and dashed blue line (i.e. Lb) are identical to those in 
Fig(3] The dashed black line (i.e. LI) is identical to that for 
SKA. 
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